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^ ■ Abstract 
(N . 

^ I We use the single particle excitation energies and the completeness rules of the 3- 

pg ■ state anti-ferromagnetic Potts chain, which have been obtained from Bethe's equation, 

T— I . to compute the modular invariant partition function. This provides a fermionic con- 

^—^ I struction for the branching functions of the Z?4 representation of Z4 parafermions which 

^vq ' complements the previous bosonic constructions. It is found that there are oscillations 

0^ ■ in some of the correlations and a new connection with the field theory of the Lee- Yang 

^ ! edge is presented. 
-)— > , 

Oh: 

^ : 1 Introduction 



X 



The theory of integrable quantum spin chains was initiated by Bethe in 1931 ||T| in his study 
of the spin 1/2 Heisenberg anti-ferromagnet. From that beginning, particularly in the last 20 
9^. ', years, an enormous number of one dimensional quantum spin systems have been discovered 
which, along with their two dimensional statistical counterparts, have the remarkable property 
that their energy eigenvalues are given by the solutions of a system of equations which have 
become known as Bethe's equations: 



j^nM+I 



sinh(Aj - iS'j) _ J^ sinh(Aj - A^ - 27) 
sinh(Aj + iS'y) j^J--^ sinh(Aj — A^ + ^7) 



Here M is the number of sites in the chain, A^ and L are related to M (typically N = M 01 
2M), and 5* and 7 are parameters which characterize the specific models. 

One of the features derived from these Bethe's equations is that in the limit M ^ 00 the 
spectrum of lowlying excitations Eex above the ground state is expressed in terms of a set of 
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single particle levels ea{P) depending on a momentum P and combined with a set of rules as 

E,,-Egs= E e-{P^), (1.2) 

a,rules 

P= Y. Pt. (1-3) 

a^rules 

and almost without exception one of the rules of combination is a "fermi" exclusion rule: 

P^^Pf if i^j. (1.4) 

The form for energy levels (1-2) and (1-3) is referred to as a quasi-particle spectrum. Further- 
more, in many of these spin chains one or more e^ vanish as P ^ 

e(P)~w|P|, (1.5) 

where v is positive and is called the speed of sound. 

Much more recently, in 1984, a powerful new formalism was invented [0 to study those 
integrable systems for which there is no mass gap and (1.5) holds. This method, known as 
conformal field theory, is more axiomatic. It deals with a continuum approximation to the 
spin chain (or two dimensional statistical system) and instead of starting from a Hamiltonian 
it starts from a symmetry principle such as the Virasoro algebra, a Kac-Moody algebra, or 
modular invariance which captures many of the essential features that lead to the integrability 
of the systems described by Bethe's equation (1.1). 

One of the main objects of computation in conformal field theory is the partition function 
which is expressed by means of modular invariance [^,^ in terms of Virasoro characters or 
(more generally) branching functions hj{q) as 

z = Y,NuM<i)H<i)- (1-6) 

Here q{q) refers to the right (left) moving excitations with 

q,q = e~^^, (1.7) 

where T is the temperature and ks is Boltzmann's constant. 

The Virasoro characters and branching functions bi{q) are solutions to the equations of 
modular transformation ||^,^,|^,^. Their construction typically starts with one or more free 
boson Fock spaces, and then excludes certain null vectors. They are typically given by explicit 
formulas with several powers of the product 

oo 

g(g) = n(l-0 (1-8) 

n=l 

in the denominator, and a power series in q in the numerator, times a fractional power which 
is usually written as g-c/24+/ifc_ Here the constant c is referred to as the central charge, and 
hfc are known as the conformal dimensions. 

The question now arises as to the relation between the solutions of Bethe's equations 
and the results of conformal field theory. In particular, one wants to compute the partition 



function (1.6) starting from Bethe's equations (1.1) (or related functional equations). Recently 



an important advance in this project was made by Kliimper and Pearce |§|,[10[! who computed 
the central charge and conformal dimensions for the Ajy+i series of the A[ models classified by 



Pasquier [jTT|. However the computation of the full character expansion (1.6) and its relation 



to the quasi-particle energy spectrum (1.2) is still lacking. 

It our purpose here to complete this project and to compute the full partition function 
(1.6) for a particular quantum spin model: The anti- ferromagnetic 3-state Potts chain. In 
particular, we will show that the partition function is constructed from the single particle levels 
of (1.2) which satisfy the fermi exclusion rule (1.4). This provides a physical interpretation 
of the model, which complements the usual computation that starts with free bosons. The 
result of conformal field theory, as obtained by Pearce [|1^], is that the partition function of 
the spin system is 

Z = J2 e-^-l^^'^ = e-^^^«/'=^^Zp;4, (1.9) 



where Cq is the ground state energy per site [|13l and Zpj4 is the D4 representation of the Z4 
parafermionic partition function of ||14||: 



Zpf4 = mq) + b'MmQ) + b'M] + mMiq) + 2bt{qK{q) + 26^(g)6^(g), (1.10) 

b^ given in one of the equivalent bosonic forms of (2.25), (2.27) and appendix B. We here 
obtain (1.10) starting from Bethe's equation (1.1) for the finite lattice and obtain new fermionic 
representations for 6^ given by (3.13) and (3.19) for 6q, fo^; andft^ in Q = and (4.16) for 6q 
and 62 in Q = =tl- 

Our method is to combine the results of [|1^, which derives the spectrum of the anti- 
ferromagnetic 3-state Potts model of the form (1.2), staring from the Bethe's equation derived 



by Albertini |T5|, with the completeness study of [|T^] and the finite size corrections of [jTO|. 
In section 2 we summarize the results of these papers which are needed here, as well as 
the conformal field theory predictions for the model. In Section 3 we compute the partition 
function in the channel Q = and in section 4 we do the same in the channel Q = 1. We 
conclude in section 5 with a discussion of the physical implication of our results in terms of 
what we call an infrared anomaly. We also discuss the oscillations which are predicted to occur 
in the correlation functions, and a connection with the field theory |T^,|T^ of the Lee- Yang 



edge [|I9| of the Ising model. 



2 Formulation 

The 3-state anti-ferromagnetic Potts chain is specified by the Hamiltonian 

2 *^ 



H 



-^ Y.{Xj + Xj + Z,Zl, + Z]Z,^,], (2.1) 



where 

Xj = / (g) / ® ■ ■ • ® X ® ■ • ■ (g) /, (2.2) 

jth 

Zj = I®I ®---® Z (g) ■ ■ ■ (g) /. (2.3) 

ith 



Here I is the 3x3 identity matrix, the elements of the 3x3 matrices X and Z are 

Xj,k = 5„fc+i (mod3), (2.4) 

Zj,k = (5j,fct^^"\ (2.5) 

uj = e2'^*/^ (2.6) 

and we impose periodic boundary conditions Zm+i = Z\. 

This spin chain is invariant under translations and under spin rotations. Thus the eigen- 
values may by classified in terms of P, the total momentum of the state, and Q, where e^'^**^/^ 
is the eigenvalue of the spin rotation operator. Here P = 2Txn/M where n is an integer 
< n < M — 1, and Q = 0, ±1. Furthermore because H is invariant under complex conju- 
gation there is a conserved C parity of ±1 in the sector Q = and the sectors Q = ±1 are 
degenerate. 

This spin chain is integrable because of its connection with the integrable 3-state Potts 
model of statistical mechanics. The eigenvalues satisfy functional equations ||15| |^ which 
are solved in terms of a Bethe's equation (1.1) [0 with 

7r/3, S = 1/4, 



iV = 2M, 



7 



and 



L = 2{M-\Q\) for Q = 0,±1. 
In terms of these A^, the eigenvalues of the transfer matrix of the statistical model are 

sinh(f ) sinh(f ) ]*^^ 
_sinh(f - A) sinh(f + A) 

the eigenvalues of the Hamiltonian (2.1) are 



A(A) 



n 



sinh(A — Afc) 



\ sinh(f| + A, 



vr 



k=l 



2M 
73' 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



and the corrsponding momentum is 

A(-i7r/12) 



JP 



^ sinh(A,. + g) 
,ysinh(A.-ft)' 



These equations have been solved to find the order one excitation energies 
are expressed in terms of three single particle excitation energies: 

e^siP) = 3(^2 cos(^ - ^) + 1}, 



IPI 



e_2.(P)=3{v^cos(y--)-l}, 



ens{P) = 3sin( 



IPI 



For P ~ 0, all three excitations are of the form (1.5) with 

V = 3/2. 
Here and in the remainder of the paper we take M to be even. 



(2.11) 
The results 

(2.12a) 
(2.12b) 
(2.12c) 

(2.13) 



1. For Q = the energies and momenta are of the form (1.2) and (1-3) 

rrta 

Ei{P!'},{P,-''},{Pn)-EGS= E T.UP^) (2-14) 



and 



where 









a=2s,-2s 


nsj-- 


=1 


p = 


a- 


>; 

=2s-2s, 


rria 
ns j=l 






pO 


= Pgs = 


M 
2" 


mod27r, 







(2.15) 



(2.16) 

m2s + m_2s is even, (2.17) 

Pj'^ ^ Pf'^'^, and Pj^'^ obey the fermi exclusion rule (1.4) and they lie in the ranges 

0<Pf <37r, (2.18a) 

< pr^' < TT, (2.18b) 

< P^' < 2n. (2.18c) 

We also note that the C parity of the ground state is 

Cgs = (-1)*'/', (2.19) 

and the C parity of an arbitrary state is 

(J / n = /I Nm„3+m_2s + (m2s+m_2s)/2 /q on") 

2. For Q = ±1 we must consider m2s + Tn^2s to be both even and odd. When m2s + ?^-2s 
is even there are two spectra of the form (2.14) and (2.15). In one P*^ = Pes and in the 
other P° = Pgs + tt. In both cases the P"" obey (2.18). When m2s + 'rn_2s is odd, there 
are again two spectra of the form (2.14) and (2.15). In each case P° = Pes- In one case 
P" satisfies (2.18), while in the other case — P" satisfies (2.18). 

The conformal field theory predictions for the model can be obtained by noting that the 
3-state Potts model is the critical D4 model in the classification of Pasquier |jTT|. The central 
charge and the conformal dimensions of the primary fields are thus obtained by specializing 
the finite size computations of the A^v+i model of Kliimper and Pearce |T0] to the case A^ = 4 



and using an orbifold construction ||2T| to find the results for D4. The general result for Ajy + 1 

(2.21) 



at the boundary of the I/II regime is that the central charge is 

2(iV-l) 



N + 2 



and the conformal dimensions are 



A^ = i(|±^-!!^ fo. H<f (2.22) 



which are the same as those of the Zn paraferinion conforinal field theory of Zamolodchikov 
and Fateev [^. Using the symmetry hi 

(2.23) 



and 



''-0 



0, 



hi 



hi 


- h^ 

— iIn 


-I 

-m 


we find for A^ = 


= 4: 




c = 


1, 








3 

4' 


hl = 


1, 


K-\ 


hi 


1 
12 



(2.24) 



where the first three conformal dimensions occur in Q = and the last two in Q = ±1. 
Moreover the modular invariant partition function is that of the D4 parafermion model |n| 
(1.10) where the branching functions 6^ can be obtained by specializing to A^ = 4 the Hecke 
indefinite form of Kac-Peterson [^ 






1(1 + 2) ni? 



Q{q) ^g 4(^^+2) 4iV 24 



X 



+ 



[ y^ y^ _ y^ y^ j /j^Ns s(s+l)/2+(«+l)n+{«+m)s/2+(Af+2)(n+s)n 

\s>On>0 s<On<0/ 

y^ y^ _ y^ y^ j {^y s{s+l)/2+{l+l)n+{l-m)s/2+{N+2){n+s) 
■,s>On>0 s<On<oJ 



for Iml < 



and using the symmetries 

bL = b' 



^m+2N 



lN-I 



(2.25) 



(2.26) 



otherwise. An alternative form for 6^ is given in |23] but for our purposes the simplest form 
is the specialization which only occurs for A^ = 4 P[ 



hl + bl- 


= /3,o/^, 6° - 6° = ^1,0/^, 6° = /3,3/2r^, 6^ = /3,2/r], 6^ 


= /3,l/^, 


(2.27) 


where 










r; = gMg(g), 




(2.28) 


and 


00 00 








fa,= E ^''^"^^^ 9a, = E (-l)V^"+^)^ 




(2.29) 




n=— 00 n=— 00 







This form has a simple origin in the Gaussian model with r = \/3/2, which we give in appendix 
A. For comparison with the expansions of subsequent sections, we list the first few terms of 
(2.27) as 

a/24L0 



gV24^o = 1 ^ ^2 ^ 2g3 + 4g4 ^ 5^6 ^ g^6 ^ ^gg^ + IQg^ + 25g^ + 37g^° + ■ 



^1/24^0 ^g3/4(l + g + 2g2 + 3^3 
-.2 I Q^3 I 



5g^ + Tq" + 12g^ 

,4 I 0^5 I 10^6 



J 



.10 



16g' + 24g» + 33g^ + 47g'^ + 



9 -2 

gi/24^0 ^ ^(1 ^ g ^ 3^^ ^ 3^^ ^ g^4 ^ g^o ^ ^3^D ^ ^7^/ ^ 27g» + 35g^ + 51g'^ + 

gV24^2 ^ ^1/3 (^^ ^ 2g + 3g2 + 5g^ + 8g^ + 13g^ + 19g^ + 28g^ + Alq^ + 58g^ + 81g^° + ■ 

gi/2452 ^ ^i/i2(^^ + g + 3g2 ^ 4^3 ^ g^4 + 11^5 ^ ^g^6 ^ 25g7 ^ 3g^8 ^ 529^ + 76g^° + 



(2.30a) 

(2.30b) 
(2.30c) 
(2.30d) 
(2.30e) 



We finally note that there is an alternative way to obtain these conformal field theory pre- 
dictions which utilizes VF4 algebra p^ and is related to the GKO construction 3 h^^\ 3 h 



i^rh 



The branching functions have been computed from this construction in terms of three dimen- 
sional sums p5|,|20,p^ which, for later comparison, we give in appendix B. 



3 Branching Functions for Q=0 

The partition function for the Hamiltonian (2.1) is, by definition, 

n 

where for M —>■ oo, 

TTC11 1 

Eos = Me. - ^ + O(^) (3.2) 



and from |TD[ c = 1. To obtain the relation with the modular invariant partition function of 
conformal field theory we must evaluate (3.1) in the limit M — >■ oo, T -^ with MT fixed. 
We intend to carry out this evaluation by making use of the quasiparticle energy spectrum 
(2.14). 

There are, however, two questions that must be addressed before we can do this. The first 
is that in order for (2.14) to specify the energy levels, the momenta P" must be discretely 
specified on the finite lattice. The second is that the evaluation leading to (2.14) is only 



correct to order one as M — i> cx) and hence in order to agree with [10] it may be necessary 
to add some term of the order of 1/M which is independent of P" but which in general will 
depend on m2s, rn_2s and rrins- These considerations are different for Q = and Q = ±1. 

We consider in this section Q = 0. We find from the previous study ||16[ of the competeness 
of the solutions of (2.1) that for given m2s, ""^-28, and rrins 



and 



Pf' takes 


oivi Iib2s -r iib^2s 1 


(3.3a) 


nins ^ \amco. 


Pr^' takes 


M ni2s + 'm_2s 

rrins values. 


(3.3b) 


P^' takes 


M — rrins — f^2s — fn-2s values. 


(3.3c) 


[ be the case if P" satisfies 








^9c TT , rrio^ + m,_29 
< P^' < Stt - —{runs + -^ + 1), 


(3.4a) 


TT , mss + m_2s , . . 


^ 9» TT , TTlos + ^71_2, 

< Pi'-' < vr - -{runs + 2 + ^^' 


(3.4b) 


ir^y^ns + "^2s + m^2s + 


1) < P/' < 27r - ^{mns + m2s + m_2s + 1), 


(3.4c) 



where the spacing between allowed values for P° is 27r/M , and m2s + Tn-2s is even. It may be 
verified that this choice of P° exactly reproduces the correct number of momenta of table 4 
of [0 for each allowed set of m2s-, Tn_2s and rrins- 



Since M ^ oo and T ^ with MT fixed, only those values of P"" where ea{P) is small 
of the order 1/M contributes to (3.1). This occurs for 

Pf ~ for a = 2s, -2s, ns, (3.5a) 

p-'^-vr, P,"^~27r ,Pf ~37r, (3.5b) 



where we may linearize ea{P) near the endpoints (3.5) as 

^a{P) ~ "W-P" for P" near zero, (3.6a) 

e-2s{P) ~ vi-n - P-^') for P^^s ^^^^ ^^ (^3 g^) 

e„,(P) ~ v{27T - P"^) for P"* near 27r, (3.6c) 

and 

e2s(P) ~^^(37r-p2") for P^' near 37r. (3.6d) 

Thus we let mf^ be the number of P" near zero and m^ be the number of P" near the end 
points (3.5b). We note that 

m^ + m^ = rria. (3.7) 

We also note that if 

m'2s + iTi^^2s is odd (3.8) 

then from (2.15) the total momentum of the state is macroscopically shifted tt from the ground 
state value Pes- These states are expected to make oscillatory contributions to the correlation 
functions. 

Consider first the case where all tti^ = (which by symmetry is identical to the case 
m^ = 0) and evaluate the partition function (3.1) using (2.14), (3.4), and (3.6) in the case 
C/Cgs = 1 where by (2.20) 

^L + ^-2s + i^L + ^-2s)/'^ is even. (3.9) 

We present in table 1 the terms from this construction up through order q^ where we see that 
they agree with the corresponding terms from the branching function g^/^'^^Q of (2.27). This 
equality has been verified to order q'^^^ and thus we conclude that this construction correctly 
gives the branching function 6°. 

To obtain a formula for 6q from the construction, let Pd{iTi, n) denote the number of distinct 
ways that the integer n can be additively partitioned into m distinct parts. Then, modifying 
the usual construction of a free fermi partition function in terms of Pd{m,n) to account for 
the momentum exclusion rule (3.4), we find 

00 00 

gi/24^0 ^ ^ ^ Pd(m„„n„,)P,(m2.,n2,)P,(m_2„n_2.)g""»+"^^+"-^^ 

rrins ,rn2s ,m-2s =0 Uns ,n2s ,n-2s =0 

m„^ , - \ ("^9c+r^._9„) , Tn9„+rri_9„ _, 

where m2s + m-2s and irins + fn-2s + (^^2s + ^^-2s)/2 are even and Pd(0, 0) = 1 by definition. 
The sums over 11^ are evaluated using 

00 „m(m+l)/2 

^P,(m,n)g" = ^-y^ , (3.11) 

n=0 \Q)m 

where we use the standard notation 

m 

(g)™= 11(1-9')' (3-12) 



and (g)o = 1 by definition. Thus we find when m2s + m-2s is even and (3.9) holds 

„m„s{m„s+l)/2 „m,2s(m2s+l)/2 „m_2s(ni_2s+l)/2 
„l/24,0 _ V- V V ^ -i ^ 

9 % - 2^ 2^ 2^ / N / X / X 

m„s=0m2a=0m_2s=0 WMns ^^"28 iyM-2s 

^(mns+m2s+m.-2s-l)„ ''"^°^2'"~^°' (mn, + (m.2s+r?i-2s)/2~l) /g j^gN 

We may now extend these considerations to the general case where both some m^^ ^ and 
some m^ 7^ 0. In this general case we note from the work of |^,|Tn| that the contribution to the 
energy from regions where (3.5a) holds and the region where (3.5b) holds are independent. 
Combining the above considerations we have in general the expression for the low lying energy 
levels in the M ^ 00 limit: 

Ee.(P/") + EejP,'-'") , (3.14) 

jf=i jf=i J 

where we define P^'" and P^'"" to satisfy 



:(mL + ^^T. '^ + 1) < ^f ^ Pt"\ (3.15a) 



^«. + < + ^'-2. + 1) < P'r% (3.15b) 

and 

^«. + ""^,"'-" + 1) < Pr\ Pr'-. (3.16a) 

^«. + ^L + ^l2s + 1) < ^^'"^ (3.16b) 

where again, the spacing between allowed values for Pj'^ is 2tt/M and ea{P) = vP with v 
given by (2.13). 

In table 2 we consider the cases mj^ = 1, rrf_2s = "^L = ^^^ m^2s = 1; "^L = "^L = 
and compute the contribution to made to Z to order q^^/^ of the terms in (3.14) that involve 
only -Pj'". From (2.17) we see that 

m2s + 'rn^-2s is odd. (3.17) 

Further, the interchange m2g <-> 'm^-2s leaves (3.14) invariant and from (2.20) gives C <-> —C. 
Thus we need only consider m2g < Tn^-^2s ^^^ fi^*^ ^^at this construction agrees with q^/'^^lf^ of 
(2.30b). 

In table 3 we consider the case rrf^^ = 1, mlg = 'rrf_2s = and compute to order q^ the 
contribution made to Z in the channel C/Cgs = 1 of the terms in (3.14) that involve only 
P/'". From (2.7) and (2.20) we find that 

TTigs + JTiigs even and m^^ + mljs + ('^L + ^-2s)/2 odd. (3.18) 

We find that this agrees with q^^'^'^b^ of (2.30c). 



These above two equalities have been verified to order q 



200 



for 


hi 


for 


hi 


for 


hi 



From these constructions we can find expressions for h\ and 64 as we did above for 6q. We 
thus find that g^^^^fe^ is given by (3.13) for a = 0, 2, 4 where 

m2s + Tn-2s 
m2s + ^«-2s is even and rUns + ^-2s H 7: is even, 

^^2s + ^^-2s is odd and m2s < fn-2s-, 

m2s + Tn-2s 
fn2s + rn_2s is even and m^^ + m_2s H is odd. (3.19) 

We may now finally construct the complete (5 = contribution to Z by using (3.14)-(3.16) 
in (3.1) and summing over all m^ and m^ subject only to the restriction (2.17) written in the 
form 

^2s + "^-2s + ^2s + ^-2s even. (3.20) 

(We note that there is no restriction corresponding to (2.20) because both channels C = ±1 
are considered in the sum.) It is easy then to see that the result consists of all the terms in 
(1.10) which involve h^ fi\ and h\ where we note: 1) that the factor of 4 in front of hl{q)hl{q) 
arises because of the symmetry under m2g ^->- ^L2s ^^^d mg^ •«-»• rn^_2si and 2) terms like 
ho{.Q)hl{q) and h\{q)hl{(i) are excluded by (3.20). 

4 Branching Functions for Q=l 

The channel Q = ±1 is more complicated that the channel (5 = because as seen in section 2 



the spectrum of excitations has 4 separate contributions. In |]I6[ these contributions are 
distinguished by the number m++ of (++) pairs of roots, and the number m__|_ of ( — (-) pairs 
of roots where there is the sum rule: 

m2s + 2m„s + ?)m_2s + fn |_ + 77i++ = M — 1. (4.1) 

We found that the three cases occurred of 

m_+-m++ = 1,0,-1, (4.2) 

and when m_j^ = m++ the spectrum is two fold degenerate. We will thus extend the consid- 
erations of section 3 by considering these three cases separately. 

4.1 m + — ni++ = — 1 

In this sector the total momentum is given by (2.15) with P° = Pas + ti" and there are M — 1 
single particle states with m„s = 1- 

We find for all three cases (4.2) from the previous study of completeness [|T6[ that 

Pj^ takes M — 1 + m_|_+ + m_2s values, (4.3a) 

P^^^ takes m_2s + "^++ — 1 values, (4.3b) 

and 

PJ"^ takes rrins + 2m++ + 2m_2s values. (4.3c) 

10 



In this present case we use r7i__|_ = 7/1++ — 1 in (4.1) to write 

2m++ = M - m2s - 2m„3 - 3m_2s, (4.4) 

and thus (4.3) reduces to 



— 1 — nins —^ values for Pf^^ (4.5a) 

2 2 ■' 

— — 1 — nins —^ values for P~^^, (4.5b) 



and 

M — nins — ""^28 — 'm-2s values for P"*, (4.5c) 

where m2s + Tn-2s is even. This will be the case if P" satisfies : 

+ 2) <P,-^-<.-^(m„. + '"" + '"-" +2). (4.6b) 

and 



M^" 


^ns 






2 




^ ^ 
M^"^ 


ns 


+ 


"^2s 


+ r/i. 


-2s 




2 





— :(m„s + m2s + "i_2s + 1) < PT < 27r - -(m„s + m2s + m_2s + 1), (4.6c) 

M ■'A 

where 

P/+1 - Pf = 27r/M. (4.7) 

It may be verified that this choice of P" exactly reproduces the momenta of table 6 of |16| 
for each allowed set of m^s-, rn^^s and m„s. Following the procedure of section 3 we compute 
in table 4 the contribution these states will make to the partition function where we use the 
linearized energies (3.6) and keep all P" near zero. 

4.2 m + — m++ = 1 

In this case P^ = Pqs and there are M — 3 single particle states with rrins = 1- Furthermore 
we find from (4.1) that 

2m++ = M - 2 - m2s - 2mns ~ 3m_2s, (4.8) 



and thus there are 



2 — nins values for Pf^ 

2 2 



-2-mns -^ ^ values for Pf, (4.9a) 

2 -^ 

- 2 - rrins values P- , (4.9b) 



and 



M - 2 - nins - m2s - m^2s values for P"*, (4.9c) 



11 



where again m2s + Tn^2s is even. This will be satisfied if P" satisfies 

TT , mas + m_2s , qn / o2s ^ o ^ / , ^a^ + m_2s , ^n ^^ m ^ 

T7(m„sH h3)<P- <37r--— (m„sH h3), (4.10a) 






^V-ns 


1 




2 




^(m„. 


+ 


m2s 


+ m_ 


-2s 




2 





m.. + :Zf^_:Zz££ + 3) < p--^s < ^ _ _(^^^ ^ ■^^^^__^^-^ ^ 3^^ ^^^^^^^ 



(m„s + m2s + m_2s + 3) < P""* < 27r - 77(^^5 + ^2^ + m_2s + 3), (4.10c) 



^(m„, + m2, + m_2. + 3) < P"" < 27r - ^, 
where (4.7) holds. Again it may be verified that this choice of P" exactly reproduces the 



momenta of table 6 of [|1^] for the allowed values of m2s, rn-2s and rUns- In table 5 we 
compute the contribution these states make to the partition function where the linearized 
energies (3.6) are used and all momenta are kept near zero. 

4.3 m + = m++ 

In this case P^ = Pes, there are 3M — 4 single particle states with 17123 = 1 and M — 4 single 
particle states with m_2s = 1. Furthermore 

2m++ = M - 1 - m2s - 2m„5 - 3m_2s, (4-11) 



where now m2s + Tn_2s is odd. In this case there is a double degeneracy and we find from |16 
that there are 



2 X (^^ - 1 - m„, - '^"^ '^"^ ) values for Pf , (4.12a) 

2 2 ■' 

■M _, m2s+m_2^ + l, , , ^_2s 



Mi ^ m2s+m_2s + l^ ______ r___ r>2s 

3 ' 

2x(— -l-m„,- '"'^^'"~^^^' ) values for P;^\ (4.12b) 

2 2 ■' 

and 

2 X (M — 1 — mns — m2s — m_2s) values for P""". (4.12c) 

Now in order to get a formula for the momentum which respects (4.12) we must consider 
two subcases. Either P" satisfy 

J^\^ns + ^ + ^)^Pj <37r-— (mn, + ^ + 3), (4.13a) 

^(m„, + ^ + l)<^j <7r-— (m„, + ^ + 3), (4.13b) 

71 71 

— :(m„^ + 17128 + m^2s + 1) < -Pr <27T - -Jji^ns + ^2^ + m_2s + 3), (4.13c) 



or P" satisfy 
- 37r + -(m„. + — ' ■■-' ' ^ + 3) < i^^ < --(m.. + "^ ' '—^ ' ^ + 1), (4.14a) 



71" ,___ , m2s + rn_2s + 1 , ON ^ r,2s ^ ^ /___ , "^2s + m_2s + 1 

J 

^ I... , "^25 + m_2s + 1 , ON ^ r.-2s ^ ^ /... , "^28 + m_2s + 1 



vr + -g(m.. + -^-- ' •I-'' ' + 3) < P,r- < --(m^^ + — ' '— ' + 1), (4.14b) 

12 



- 27r + ^{rrins + m2s + m_2s + 3) < P"' < -^{rUns + m2s + m^2s + !)• (4.14c) 

M ■' IVl 

However, now instead of the total momentum P being given in terms of P° by (2.15) we 
must introduce a shift of order 1/M (which is permissible because (2.15) is only derived to 
order one) and write that when (4.15) holds 

P = P'+ E !:^° + ^( ""^"-""' ), (4.15a) 



a=2s,-2s,nsj=l " ^'^ 



and when (4.16) holds 



a=2s-2s,ns j=l ^"^ ^ 

It may be verified that the momenta computed from these rules agrees with the momenta of 
table 6 of [0 for the allowed values of m2s, m_2s and rrins- 

Corresponding to this momentum shift there is an energy shift as well. Thus in table 6 we 
compute the contribution to the partition function of those states obtained from (4.13) with 
P" near zero using the linearized energies (3.6) and subtracting the shift ■iL( ™23+r»-2s- y ^^^ 
note that the macroscopic momentum of these states are near Pes- In table 7 we compute the 
contribution to the partition function from those states obtained from (4.14) with Pj^ near 
— 37r, Pj^^* near — vr and P^'^ near — 27r where we add a shift ■f^i "^'^" "^~'^''~ ) to the energy. 
The macroscopic momentum of these states is Pes + tt. In both cases, due to the symmetry 
under 2s ^-> —2s, only the states with m2s < fn-2s are shown. 

4.4 Branching functions 

We may now obtain the formulas for the two branching functions oiQ = 1, namely 63 and b"^ by 
combining together the results of the three preceding subsections with the same macroscopic 
momentum. 

Consider first the case where the macroscopic momentum is Pes- This is obtained from 

the m ^ — m++ = 1 states of table 5 and the m__|_ = m_|_+ states of table 6. In table 8 we 

compute the sum of these two contributions and see that it is identical with the corresponding 
terms in q-^/^q^/^%1 of (2.30d). 

The other case has the macroscopic momentum of Pes + ^r .This is obtained from the 

m ^ — m^^ = — 1 states of table 4 and the m ^. = m^^ states of table 7. In table 9 we 

compute the sum of these two contributions and see that it is identical with the corresponding 
terms in q'^/^^q^/^^l of (2.30e). 

As in Q = these identities have been verified to order q"^^^. 

We may now use the above construction to obtain formulas for the branching functions 6q 
and 62 by following exactly the same procedure used in section 3. Thus we find 
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„ „ „ mnsjmns+l) "I2s('"2s + 1) "^-23 ('"-2a +1) 

°° °° °" g 2 g 2 g 2 



m„s=0m2s=0 ™-2s=0 K'djrrins \H)m2s \H)m-2s 

m,2s+m-2s even 

m^o / -. \ ('"9«+'" — 9c) / ^T19„+m_9„ _, 

^ ^(mna+m2s+m_2a+l) ^ 2 ^('""3 + ^^^ 2£+i) 

„ „ „ mn3(mn3+l) ™2s('^2s+l) '"-2s ("^-2s +1) 

g 2 g 2 g 2 



+ E E E ^^ 

m„s=0m2s=0 m-2s=0 [Q)m„3 [Q)m2s [Q)m^2 

m2s+m-2s odd 



'"ns 



, ^, (m9„+Tn_9„), m9a+m_9„, _ , 

(mns+m2s+m_2s-l)„ ^° 9 2£i(m„^ + ^^^ i±) i/^ 



and 



^ ^^(m„s+m2s+m-2s-ljg^ ^(m„s + ^^:^^ ^)„-i/4 (•4.1ga 



^ „ ^ mnsi.mns + i) '"2s('"2s + l) '"-2s ('"-23 +1) 

°° °° °° g 2 g 2 g 2 



,-"","-^1 = E E E ,,) ,,) ,,) 

m„s=0m2s=0 m_2s=o \H)mns \H)m2s WM- 

m.2s+m-2s even 

m„. / , I -,^ ('"2s+'"-2s) / I '"2s+'"-23N 
^ „— f^(m„s+m2s+m-2s-l)„ 9 (m^sH 9 ) 



^ „ „ '"n3('"ns + l) '"2s('"2s+l) '"-2s ('"-2s +1) 

m„s=0m2s=0 '"-2s=0 \^)m„s lgjm.2s \^)m-2a 

m2s+m-2s odd 

^^^(m„s+m2s+m-2s+l)g ' 2 ("ins+ 9 ^° +1)^1/4^ (4.16b) 

4.5 Partition function 

Finally we need to compute the complete partition function in the Q = 1 channel. Now in 
each of the four cases considered above there are both left and right excitations m^ and m^. 
Denote the four sums obtained above for all mj^ = by S*" for m_+ — m^^ = — l,by S~^ 
for m_+ — m++ = 1, by S*" for m_+ = m++ with AP = 0, and by S'^ for m_+ = m^^ 
with AP = 71. To consider the general case with both m^ and m^ both nonzero, we follow 
the procedure of Section 3 and consider momentum restrictions for the right and left movers 
separately. 

Consider first m_+ — m++ = —1. Then because of the restriction that m2s + m_2s = 
m2s + Trf^2s + "^-2s + "^-2s ^ust be even we see that m2s + ^!l2s ^^^ '^L + "^-2s must be even 
or odd together. If both terms are even the contribution made to the partition function is 
S^{q)S~{q) and if both terms are odd and the momenta are shifted properly the contribution 
is S^{q)S^[q). 

The contribution from m |_ — m++ = 1 is similar. Again nf^g + fTf-2s ^"^^ ^s + ^^-2s uiust 

be even or odd together. The even terms contribute >S'"'"(g)5'"'"(g) and the odd terms contribute 
S\q)S\q). 

Finally there are the terms from m++ = m^j^. Here m2g + Trf-2s + ^^L + ^-2s is odd, thus 
777.25 + "^-2s is even (odd) and 7772^ + Tn^^2s is odd (even). The ends where 7772s + ^-2s ^'^^ odd 
produce the terms S^ and S'^ and if the momentum shift is properly accounted for the ends 
with 77725 + ^-2s even produce S^ and S~ . Thus the sector 777++ = 7r7_+ produces the cross 
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terms like S'-' (q) S'^ (q) and S'^{q)S (g) and hence the desired form bl{q)bl{q) + bl{q)bl{q) in 
(1.10) is obtained. 

5 Discussion 

There are many mathematical points of the foregoing computations which remain to be clar- 
ified such as: 1) a direct proof of the equivalence of the forms (3. 13), (3. 19) and (4.16) with 
(2.25), (2.27) and the forms of appendix B; 2) the obtaining of our results directly from the 



functional equations [|T5|,0 without recourse to the study of the completeness rules |jT^. How- 
ever the major feature of the results of sections 3 and 4 is that they directly relate the concept 
of branching function with that of quasi-particle. This provides insight into the physics of the 
model which we will discuss here in detail. 

5.1 Infrared Anomaly 

All the eigenvalue spectra computed here have, up to a possible constant, the additive quasi- 
particle form (1.2) of 

Ee.-EGS= E ea{Pj), (5.1) 

a,rules 

where there are two important features of the rules which govern the combination of energy 
levels. The first is the fermi exclusion property: 

P^'-^P- for t^J. (5.2) 

This rule in conjunction with the quasi-particle form (5.1) is often used to say that the quasi- 
particles are fermions. 



If the momenta P" were such that 



2vrj 
M 



P^ = ^ (5-3) 



with j an integer (or possibly half integer) these quasi-particle energies would indeed be 
identical with those of a genuine free Fermi gas. However the momentum set out of which 
the Pj are chosen is not (5.3). Instead the momenta are subject to the restrictions (3.16), 
(4.6), (4.13) and (4.14). These all share the feature that there is a depletion of the number 
of allowed momenta near the values of P where ea{P) = 0, which depends on the number of 
quasi-particles in the state. This is an intrinsic many body effect which cannot be modeled 
by an effective change in some one or two body property. Because this occurs for e(P) ~ 
it seems appropriate to refer to this intrinsic many body effect as an infrared anomaly. This 
mechanism is typical of all conformal field theories that are different from free fermions or free 
bosons and is what can lead to central charges being different from integer or half integer. 

This manner of characterizing the phenomena underlying the behavior of systems with 
nontrivial branching functions has, at least on the face of it, nothing to do with integrability, 
Virasoro algebra, modular invariance, or any other symmetry algebra. In this respect it shares 



the feature of generality with Haldane's |]28| definition of fractional statistics. However the 
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definition of |^ differs from that of tlie infrared anomaly by relying on tlie finiteness of tlie 
Hilbert space. Tliis can only be achieved by imposing an ultraviolet as well as an infrared 
cutoff on the problem whereas, by the very name, an infrared anomaly will exist without 



an ultraviolet cutoff. The essential feature of [^, however, is the abandoning of a second 
quantized description of the excitations in the system and this is certainly a key property of 
the effects described above. 

5.2 Specific Heat 



One of the powerful results of conformal field theory is the prediction |29|,Q that the specific 
heat C, for a system with periodic boundary condisions, is given in terms of the central charge 
and the velocity of sound as T ~ as 

C ^ ^T. (5.4) 

The original derivations are based on conformal invariance. We give here an argument from 
the point of view of the infrared anomaly. 

By definition the bulk free energy per site is obtained from the partition function 

Z = Tre"^ (5.5) 

as 

/ = -ksT lim —InZ, (5.6) 

where 

T > is fixed and M ^ oo. (5.7) 

Then the specific heat is 



C = -ksT^,. (5. 



,5V 

To obtain the leading behavior as T ^ of the specific heat it is sufficient restrict attention 
to the lowlying (order 1 or smaller) excitations of H over the ground state. These obey the 
quasi particle form (5.1). Hence specific heats are commonly evaluated with formulae involving 
single particle levels e(P). In this discussion the boundary conditions are involved only in a 
tacit fashion. 

This argument, however, is not complete, as is apparent in the observation that any energy 
level with e(P) > as M — i> cx) will contribute only a term exponentially small in T to the 
specific heat. Thus the order one excitations do not contribute to the linear term (5.4). Instead 
it is the levels which have the property that lim^^oo e(-P) = that contribute to the leading 
behavior. 

The partition function computed in conformal field theory is in the limit 

M ^ oo, T ^ with MT fixed. (5.9) 

This is not the same as the limit (5.7) which defines the specific heat. However if no additional 
length scale appears in the system it is expected that the behavior of the specific heat computed 
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using the prescription (5.9) where q = e '^"'bt is fixed will agree when q ^ 1 with the T ^ 
behavior computed using the prescription (5.7). 

The g — » 1 behavior of Z can be computed directly from the expression of the branching 
functions as infinite series in q. Since these series are a direct consequence of the infrared 
anomaly the infiuence of the many body effects which determine the momentum selection 
rules is apparent. These selection rules can only be seen by imposing an explicit infrared 
cutoff on the problem. Indeed it is known [^ that free boundary conditions and periodic 
boundary conditions give different amplitudes of the linear in T term of the specific heat. 

The above argument has relied only on a one length scale scaling argument and has not 
explicitly used modular invariance. Modular invariance in this context is a mathematical 
property of the series expressions for the branching functions which says that if we define r 
as g = e^'^" then for a set of branching functions hk{T) 

bk{-l/r) = J2MkMr), (5.10) 

where the functions Mk^i are independent of r. Thus the q —>■ 1 behavior of the branching 
functions is given in terms of the g — > behavior of the branching functions. This behavior 
is always g'M+^fe where hk is the conformal dimension and c is determined from the finite 
size correction to the ground state energy of (3.2). The right hand side of (5.10) will thus 
be dominated by the term with the smallest h^. In particular if the smallest h^ is zero the 
formula (5.4) results. 

This argument has now related the specific heat to the finite size correction to the ground 
state energy and here it is apparent by definition that boundary conditions are important. 
Nevertheless it is useful to keep in mind the fact that the specific heat would not be affected 
if we changed the finite size corrections provided that the many body effect of the infrared 
anomaly was not changed. In that case the modular invariance would be destroyed but the 
specific heat would be unchanged. This is a demonstration that the concept of infrared 
anomaly and modular invariance are logically distinct. 

5.3 Oscillations 

One of the striking features of sections 3 and 4 is the fact that the branching functions 63 
(with conformal dimension 3/4) and 63 (with conformal dimension 1/12) are obtained with a 
macroscopic (order 1) momentum shift from the ground state of AP = it. This results from 
the feature found in [|T3| that the energies e2s{P) (e_2s(-P)) vanish at 37r and vr. These macro- 
scopic momentum shifts are expected to give rise to oscillatory contributions to the correlation 
functions of the primary operators of h\ and h\. On the lattice the oscillatory term should be 
(—1)^ where A^ is the separation of the operators. These microscopic oscillations are, perhaps, 
unusual for conformal field theories but are in fact expected if we make the observation that 
the anti-ferromagnetic 3-state Potts model lies deep inside the incommensurate phase of the 



chiral Potts model |15, 31 1 which is characterized by oscillatory correlations ||32|, 
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5.4 Lee- Yang Edge 

An interesting property of the branching functions (3.13) and (4.16) obtains if we consider 
only the terms where m2s = iTi^2s = where the sums reduce to the two sums 

So= Y. f—- (5.11a) 



and 



These are the famous sums of Rogers- Ramanuj an P^ , pB| . They become modular functions if 
we multiply by powers of q and consider 



oo ^mns(ni„s+l) 

^1= E '-?-. ■ (5-llb) 

m„s=0 \H)m„s 



ci,3{r) = q ^0 So and ci,i(r) = geo^^, (5.12) 



which, in fact, are shown in [Q to be the two characters of the nonunitary minimal model 
with p = 5 and p' = 2 

-1/24 oo 
^ rg\ — 1 y^ I (2npp'+rp-sp')^/4pp' _ (2rapp'+rp+sp')^/4pp'| /g j^gN 

Q[Q) n=-oo 



This model has been identified [17, 18| with the field theory that describes the behavior at 



the Lee- Yang edge of the Ising model [|T9[. Thus there is a sense in which we may say that 
the field theory for the Lee- Yang edge is obtained by adding a perturbation to the 3-state 
anti-ferromagnetic Potts chain which makes the ±2s excitations massive without affecting the 
ns excitations. 

It is also interesting to note that the g — ;> 1 behavior of the sums (5.11) can be calculated 
directly without recourse to the modular transformation (5.10) by a method that makes con- 
tact with dilogarithms. As an example we consider explicitly 5*0 which we write in the form 
obtained directly from (3.10) 

-^0 = E E P'ii^^ n)g"gT(— 1) (5.14) 

m=0 n=0 

(where the subscript ns has been dropped for simplicity). To study the limit g ^ 1 we first 
use the integral representation 

j:nm,n)q- = —f-^l[il + zq^) (5.15) 

in the sum (5.14). We note that (5.15) vanishes if ??7, < 0. Thus we extend the lower limit of 
the sum over m from zero to — cxd and interchange the summation and integration to obtain 



5o = ^/-n(l + W) E g-^"-^)^-. (5.16) 



dz °° 

1=1 m=—oo 



The sum over in is expressed in terms of the Jacobi theta function 

oo 



m=—oo 



as 



with 
and hence 



g gf(™-l)2— = g-i^l02(t^,g^), (5.18) 

m=— oo 

e^^^ = z-|, (5.19) 

^0 = ^ / ^ n(l + zq')q-hk,{v, q^). (5.20) 

Then using the product representation 

oo 
^2(t;,gl)=gi(^l+^-|) J](l-g")(l + g"^)(l + g"^-l), (5.21) 



n=l 

we find 



S'o = exp|f;in(l-g")|^^^(^ + l)exp|f;{21n(l + ^g") + ln(l + z-ig")}|. (5.22) 
We may now study behavior of S'o as g ^ 1 by replacing the sums in (5.22) by integrals. 



Thus using the definition q = e ^-^''bt and setting 



- = 1^ (5-23) 



we find 



^ \MkBT /■» , , ,^ A 

So ~ exp<^— / rfa;lnl-e^U 

[ 271V Jo J 

1 r rl7 ( MIctdT r°° ^ 

X _ i _(^ + 1) exp i -^^^^^ / dx{2 ln(l + ze-"") + ln(l + z-^e-"")} K5.24) 
27ri J z [ 2'Kv Jo J 

The integral over z may now be evaluated by steepest descents. The steepest decents point 
occurs at the values of z that satisfy 

ln(l + 2)2 = ln(l + z-^), (5.25) 

and thus either 2; = — 1 or 

l^z = z'\ (5.26) 

and hence we find that the steepest descents point is 

. = ^. (5^27) 
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Thus we have 



So ~ exp <^ -—^- / dx{\n{l - e"") + 2 ln(l - 
[^ znv Jo 

which, if we define 

t = e-^ 

and recall the definition |3^ of the dilogarithm 



V5- 
7. 


1 
— e 


'") + ln(l + 


V5 + 1 
2 ' 


(5.28) 


















(5.29) 



Li2(z) = - r dt^^!^^—-^, (5.30) 

JO t 



may be rewritten as 



we obtain the result 



5*0 ~ exp <^ ^— — {Li2(l) + 2Li2( — ^ — ) + Li2( ^ — )} \ . (5.31) 

Then if we note the special values of the dilogarithm 

Li2(l) = ^, (5.32a) 

o 

L,(-i±i^,^-^-.„^(^). (6.320 

MkeTn 

So ~ exp -—- . (5.33) 

3Uf 

The identical behavior is obtained for 5*1. Thus with Z = So{q)So{q) + Si{q)Si{q), we obtain 
from (5.4)-(5.8) an (effective) central charge of 2/5 which agrees with ||18|| . 
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Appendix A. Gaussian construction of the branching func- 
tions 

N The form of the branching function (2.27) may be simply obtained if we note that the 
Z4 parafermions with the diagonal {A^) modular invariant partition function is known to be 



the r = \/3/2 point on the orbifold line of c = 1 conformal field theories (see e.g. |3^ 
However, we are interested in the nondiagonal (-D4) partition function which contains the 
current operator of dimension (1,0) in the spectrum. This model must therefore lie on the 



c = 1 Gaussian line at the compactification radius r = \/3/2. Indeed, consider the general 
Gaussian model partition function 

-| 00 

where 

A™,„(r) = -(^ + nrf and A^,„(r) = ^{^ - nrf (A.2) 



and set r = \/3/2 to obtain 



zJ^) = —l— V qH^^rf(^^)\ (A.3) 



m,n=—oo 



Rewriting the sum as 



^V2^ viqHq)to ^t-^ 

m+3n=bmod 6 



we see that 

^^§-l^,t.hMhM (A.5, 

with /m,n defined by (2.29). Then noting the symmetry 

/3,1 = /3,5 and /3,2 = /3,4 (A.6) 

we obtain precisely the Zpj4 of (1.10) with the expressions of (2.27) for 6q + 6°, 62) ^O' ^"^^ ^2- 

Appendix B: Branching functions of (A3 )i x (A3 ^)i/(A3 ')2 

The branching functions (2.27) can also be expressed as a three dimensional sum in the 
following way ||2^,^^. Let a^, z = 1, 2, 3 be the simple roots of A-^. Then 

(ai,aj) = Cjj, (B.l) 
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where c is the Cartan matrix of A-^ 







-I 



-1 



fB.2) 



Let Aj be the fundamental weights of A^, so that \j = Y.k'^jl^k- The dominant weights a*^*^) 
of level k are defined to be: 



Ak) 



Y^ aiXi, tti G Z, aj > 0, ^ai < k. 



(B.3) 



i=l 



i=l 



Let r = a*^^^ + p, s = a*^^) + p, where p = J2i^i- Then the following branching functions are 
identical with (2.27): 

&r,s = - E E det(ti;)g 60 (B.4) 



where W is the Weyl group of A^, generated by the three simple reflections: 

(Ti{P) = P - {ai, P)ai, i = 1,2,3 



(B.5) 



This group has 24 elements made up of powers of the simple reflections, and det(w) = ±1, 
depending on whether the minimum number of simple reflections making up the element w is 
even or odd. Q is the root lattice of A3, i.e. k = J2imiai, and the sum X^keQ i^ thus a sum 
over rrii, i = 1,2, 3, from —00 to 00. 

Equation (B.4) gives 7 unique branching functions, but only 5 appear in our model: these 



can be obtained by setting r = p and choosing the following a' 



(2) 



s-p: 



,(2) 



(0, b, 

2A2, K 

2Ai, b, 

A2, h 

Ai + A3, 61 



/)° h 

bi K 

hi K 

bl, h, 

bl, h, 





1 

3/4 

1/12 

1/3. 



(B.6) 



Here, the conformal dimension h^^s is defined to be: 

I — 6r + 5sP 



hy 



I 

12 



60 



(B.7) 



The three dimensional sum (B.4) is to be compared with the expressions (3.13) and (4.16) 
of the text. We note that the sum (B.4) has a power q^^^"^] and thus the powers of q grow 
much more rapidly as a function of nij than do (3.13) and (4.16). Note also that when the 
sum over W is performed (B.4) contains 24 triple sums whereas (3.13) has only one and (4.16) 
has two. 
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Table 1: The terms through order q^ in the construction of b^ from the rules of section 3. The 



ffirrins + m2s + m_2s + 1) and P, 



±2s 
ns 



^im. 



+ 



m23+m_2s 



minimum momenta P;;^^„ - j^^.^ns . ■■■^zs . ■■—zs . ^j ^^^^ ^ ns Mv-ns . 2 

obtained from (3.15). The terms in g^/^^feg are obtained from (2.30a). Here 11123 
m2s+rn-23 ^^^ eveu. The macroscopic momentum shift is AP = 0. 



+ 1) are 
m_2s and 



m. 



+ m_2s + 



order 


"^L 


"^L 


^-2s 


pns 

min 


p±'2s 
min 


|pns^ p2.^ p-2s| (^j^j^g of ^) 


states 


^1/24^0 


q^ 











— 


— 


{0,0,0} 


1 


1 


q' 





1 


1 


— 


2n/M 


{0,2,2} 


1 


1 


q' 





1 


1 


— 


27r/M 


{0,4, 2}, {0,2, 4} 


2 


2 


q' 



2 


1 



1 




37r/M 


27r/M 


{0,6, 2}, {0,4, 4}, {0,2, 6} 
{3 + 5,0,0} 


3 

1 


4 


q' 



2 


1 



1 




37r/M 


27r/M 


{0,8,2},{0,6,4},{0,4,6}, 

{0,2,8} 

{3 + 7,0,0} 


4 
1 


5 


q' 




2 
1 
1 


1 


2 



1 



2 


37r/M 
47r/M 
47r/M 


27r/M 

37r/M 
37r/M 


{0,10,2},{0,8,4},{0,6,6}, 

{0,4, 8}, {0,2, 10} 

{3 + 9, 0,0}, {5 + 7, 0,0} 

{4,3 + 5,0} 

{4,0,3 + 5} 


5 
2 
1 
1 


9 


q' 




2 
1 
1 


1 


2 



1 



2 


37r/M 
Ati/M 
A-k/M 


2ti/M 

371/M 
Sir/M 


{0,12,2},{0,10,4},{0,8,6}, 
{0,6, 8}, {0,4, 10}, {0,2, 12} 
{3 + 11, 0,0}, {5 + 9, 0,0} 
{6, 3 + 5,0}, {4, 3 + 7} 
{6, 0,3 + 5}, {4, 0,3 + 7} 


6 
2 
2 
2 


12 


g« 




2 
1 
1 



1 


2 


2 


1 



2 

2 


3n/M 
A-k/M 
A-k/M 


27r/M 

3n/M 
Stt/M 
37r/M 


{0,14, 2}, {0,12, 4}, {0,10, 6} 

{0,8, 8}, {0,6, 10}, {0,4, 12} 

{0,2,14} 

{3 + 13, 0,0}, {5 + 11, 0,0}, 

{7 + 9,0,0} 

{8, 3 + 5,0}, {6, 3 + 7,0}, 

{4, 3 + 9,0}, {4, 5 + 7,0} 

{8, 0,3 + 5}, {6, 0,3 + 7} 

{4, 0,3 + 9}, {4,, 0,5 + 7} 

{0,3 + 5,3 + 5} 


7 
3 

4 

4 
1 


19 



23 



Table 2: The terms through order g^^/^ in the construction of 62 from the rules of section 3. 



The minimum momenta P, 



ns 
min 



M 



[rrir, 



+ m2s + m_2s + 1) and P^f^ 



M 



K^^ns + 



m2s+m_2s 



+11 



are obtained from (3.15). The terms in g^'^^62 ^-re obtained from (2.30b). Here m2s 
odd and m2s < itl^2s- The macroscopic momentum shift is AP = ir. 



m-2s IS 



order 


"^L 


mi. 


m^-2s 


pns 

min 


p±'ls 
min 


|pns^ p2s^ p-2s| (^j^j^g of ^) 


states 


^1/24^0 


^3/4 








1 


— 


37r/2M 


{0,0,3/2} 


1 


1 


^7/4 








1 


— 


37r/2M 


{0,0,7/2} 


1 


1 


^11/4 




1 






1 
1 


3n/M 


37r/2M 
57r/2M 


{0,0,11/2} 
{3,0,5/2} 


1 
1 


2 


gib/4 




1 






1 
1 


3ti/M 


37r/2M 
57r/2M 


{0,0,15/2} 

{3, 0,9/2}, {5, 0,5/2} 


1 
2 


3 


^19/4 




1 







1 


1 
1 

2 


3ti/M 


37r/2M 
57r/2M 
57r/2M 


{0,0,19/2} 

{3, 0,13/2}, {5, 0,9/2}, {7, 0,5/2} 

{0,5/2,5/2 + 9/2} 


1 
3 

1 


5 


^23/4 




1 







1 


1 
1 

2 


371/M 


37r/2M 
57r/2M 

57r/2M 


{0,0,23/2} 

{3, 0,17/2}, {5, 0,13/2} 

{7, 0,9/2}, {9, 0,5/2} 

{0, 9/2, 5/2 + 9/2}, {0, 5/2, 5/2 + 13/2} 


1 

4 
2 


7 


^27/4 




1 



2 






1 





1 
1 

2 

1 
3 


Stt/M 
in/M 


37r/2M 
57r/2M 

57r/2M 

7tt/2M 
5n/2M 


{0,0,27/2} 

{3, 0,21/2}, {5, 0,17/2}, {7, 0,13/2} 

{9, 0,9/2}, {11, 0,5/2} 

{0, 13/2, 5/2 + 9/2}, {0, 9/2, 5/2 + 13/2} 

{0, 5/2, 5/2 + 17/2}, {0, 5/2, 9/2 + 13/2} 

{4 + 6,0,7/2} 

{0,0,5/2 + 9/2 + 13/2} 


1 

5 

4 
1 
1 


12 


^31/4 




1 



2 






1 





1 
1 

2 

1 
3 


371/M 
Ati/M 


3n/2M 
5tt/2M 

5n/2M 

771/2M 
57r/2M 


{0,0,31/2} 

{3, 0,25/2}, {5, 0,21/2}, {7, 0,17/2} 
{9, 0, 13/2}, {11,0, 9/2}, {13, 0, 5/2} 
{0, 5/2, 5/2 + 21/2}, {0, 5/2, 9/2 + 17/2} 
{0, 9/2, 5/2 + 17/2}, {0, 9/2, 9/2 + 13/2} 
{0, 13/2, 5/2 + 13/2}, {0, 17/2, 5/2 + 9/2} 
{4 + 8, 0,7/2}, {4 + 6, 0,11/2} 
{0,0,5/2 + 9/2 + 17/2} 


1 
6 

6 
2 
1 


16 



24 



Table 3: The terms through order q^ in the construction of 64 from the rules of section 3. The 



minimum momenta P" 



igim. 



+ 17123 + "^-2s + 1) and P^j, 



M^^r. 



+ 



m2s+m-2s 



are obtained from (3.15). The terms in g^/^^64 are obtained from (2.30c). Here m2s + Tn_2s is 



even and rUns + "^-2s + 



m2a+m^2s 



is odd. The macroscopic momentum shift is AP = 0. 



order 


"^ns 


"^2s 




pns 

min 


-'■min 


{P"^ P2^ p-2^} (units of ff) 


states 


q"''% 


q' 


1 








2n/M 


— 


{2,0,0} 


1 


1 


q' 


1 








27r/M 


— 


{4,0,0} 


1 


1 


q' 


1 





2 





2 


27r/M 


2ti/M 
27r/M 


{6,0,0} 
{0,2 + 4,0} 
{0,0,2 + 4} 


1 
1 
1 


3 


q' 


1 






2 





2 


2n/M 


27r/M 
2ti/M 


{8,0,0} 
{0,2 + 6,0} 
{0,0,2 + 6} 


1 
1 
1 


3 


q' 


1 




1 



2 

1 




2 
1 


27r/M 
A-k/M 


2ti/M 
2ti/M 
37r/M 


{10,0,0} 

{0,2 + 8,0}, {0,4 + 6,0} 

{0,0,2 + 8}, {0,0,4 + 6} 

{4,3,3} 


1 
2 
2 
1 


6 


q' 


1 




1 



2 

1 




2 
1 


27r/M 
47r/M 


2tx/M 
2n/M 
37r/M 


{12,0,0} 

{0,2 + 10,0}, {0,4 + 8,0} 
{0,0, 2 +10}, {0,0, 4 + 8} 
{4, 5, 3}, {4, 3, 5}, {6, 3, 3} 


1 
2 
2 
3 


8 


q' 


1 




1 



2 



1 





2 

1 


27r/M 
47r/M 


27r/M 
27r/M 
37r/M 


{14,0,0} 

{0,2 + 12,0}, {0,4+10,0}, 

{0,6 + 8,0} 

{0,0,2+ 12}, {0,0,4+ 10}, 

{0,0,6 + 8} 

{4, 7, 3}, {4, 5, 5}, {4, 3, 7} 

{6, 3, 5}, {6, 5, 3}, {8, 3, 3} 


1 
3 
3 
6 


13 


g« 


1 




1 



2 



1 





2 

1 


2tx/M 
AixjM 


2tx/M 
2tx/M 
37r/M 


{16,0,0} 

{0,2 + 14,0}, {0,4 +12,0}, 

{0,6 + 10,0} 

{0,0, 2 +14}, {0,0, 4 +12}, 

{0,0,6+10} 

{4, 9, 3}, {4, 7, 5}, {4, 5, 7} 

{4, 3, 9}, {6, 7, 3}, {6, 5, 5} 

{6, 3, 7}, {8, 5, 3}, {8, 3, 5} 

{10,3,3} 


1 
3 
3 

10 


17 



25 



Table 4: The terms through order q^ in the sector m_+ — m_|_^ 
the rules of section 4.1. The minimum momenta P, 



[rrir 



+ 



m2s+m-2s 



p±2s ^ 

-'■min j^V'ns ' 2 

macroscopic momentum shift is AP 



M 



.rrir 



= —1 constructed from 

+ m2s + ■m_2s + 1) and 

+ 2) are obtained from (4.6). Here m2s + itl_2s is even and the 

TT. 



ns 
min 



order 


^L 


"^L 


"^-2s 


pns 

min 


p±-2s 

^ min 


|pns^ p2s^ p-2s| (^^j^i^g Qf ^) 


states 


total 


g" 











— 


— 


{0,0,0} 




1 


q' 


1 








2n/M 


— 


{2,0,0} 




1 


q' 


1 








2n/M 


— 


{4,0,0} 




1 


<l' 


1 






1 




1 


27i/M 


37r/M 


{6,0,0} 
{0,S,S} 




2 


q' 


1 


2 






1 


2 





1 



2 


2tx/M 
Stt/M 


Stt/M 

Stt/M 
Stt/M 


{8,0,0} 

{0,S,5},{0,5,S} 
{S + 5,0,0} 
{0,S + 5,0} 
{0,0,S + 5} 


2 


6 


q' 


1 

2 






1 


2 




1 


2 


2tx/M 
Stt/M 


Stt/M 

Stt/M 
Stt/M 


{10,0,0} 

{0,S,7},{0,5,5},{0,7,S} 
{3 + 7,0,0} 
{0,S + 7,0} 
{0,0,S + 7} 


S 


7 


q' 


1 


2 


1 



1 


2 

1 



1 



2 
1 


27r/M 
37r/M 
Atx/M 


Atx/M 


{12,0,0} 

{0,S,9},{0,5,7}, 

{0,7,5},{0,9,S} 

{S + 9, 0,0}, {5 + 7, 0,0} 

{0,S + 9,0},{0,5 + 7,0} 

{0,0,S + 7},{0,0,5 + 7} 

{4,4,4} 


4 
2 
2 
2 
1 


12 


q' 


1 


2 


1 
1 
1 



1 


2 

1 
2 




1 



2 
1 

2 


27r/M 

37r/M 

Atx/M 
Ati/M 
Att/M 


Stt/M 

Stt/M 
Stt/M 

47T/M 

4tt/M 
4tt/M 


{14,0,0} 

{0,S, 11}, {0,5, 9}, {0,7, 7}, 

{0,9,5},{0,11,S} 

{S + 11, 0,0}, {5 + 9, 0,0} 

{0,S + 11,0}, {0,5 + 9,0} 

{0,0,S + ll},{0,0,5 + 9} 

{4, 4, 6}, {4, 6, 4}, {6, 4, 4} 

{4,4 + 6,0} 

{4,0,4 + 6} 


1 

5 
2 
2 
2 
S 
1 
1 


17 


q' 


1 


2 


1 

1 
1 



1 


2 

1 

2 




1 



2 
1 


2 


2ti/M 
3n/M 

47r/M 

47r/M 
Att/M 


Stt/M 

Stt/M 
Stt/M 
4tt/M 

4tt/M 
4tt/M 


{16,0,0} 

{0,S, IS}, {0,5, 11}, {0,7, 9}, 

{0, 9, 7}, {0,11, 5}, {0, IS, S} 

{S + IS, 0, 0}, {5 + 11, 0, 0}, {7 + 9, 0, 0} 

{0, S + IS, 0}, {0, 5 + 11, 0}, {0, 7 + 9, 0} 

{0, 0, S + IS}, {0, 0, 5 + 11}, {0, 0, 7 + 9} 

{4, 4, 8}, {4, 8, 4}, {8, 4, 4}, 

14 6, 6}, {6, 4, 6}, {6, 6, 4} 

|l,4 + 8,0},{6,4 + 6,0} 

{4, 0,4 + 8}, {6, 0,4 + 6} 


1 

6 
S 
S 
S 

6 
2 
2 


26 



Table 5: The terms through order 

rules of section 4.2. The minimum momenta P, 



q^ in the sector m_+ 



m 



m2s+m-2s 



M 



irrir, 



M^'^ns+ '"""'^2" ^^ +3) are obtained from (4.10). Herem2s + ^W- 



= 1 constructed from the 
17123 + "^-2s + 3) and P^f^ = 
2s is even and the macroscopic 



momentum shi 


ft is AP = 0. 












order 


^ns 


mi. 


mUs 


pns 

rain 


p±2s 


|pns^ p2s^ p-2s| (^j^i^g Qf ^) 


states 


total 


q' 











— 


— 


{0,0,0} 




1 


q' 


1 








An/M 


— 


{4,0,0} 




1 


q' 


1 








Att/M 


— 


{6,0,0} 




1 


q' 


1 








Att/M 


— 


{8,0,0} 











1 


1 


— 


4tt/M 


{0,4,4} 




2 


q' 


1 








An/M 


— 


{10,0,0} 











1 


1 


— 


47T/M 


{0,4, 6}, {0,6, 4} 











2 





— 


4tt/M 


{0,4 + 6,0} 














2 


— 


4tt/M 


{0,0,4 + 6} 




S 


q' 


1 








47r/M 


— 


{12,0,0} 











1 


1 


— 


4tt/M 


{0,4, 8}, {0,6, 6}, {0,8, 4} 











2 





— 


4tt/M 


{0,4 + 8,0} 














2 


— 


47T/M 


{0,0,4 + 8} 








2 








Stt/M 


— 


{5 + 7,0,0} 




7 


q' 


1 








An/M 


— 


{14,0,0} 











1 


1 


— 


4tt/M 


{0,4, 10}, {0,6, 8}, 
{0,8, 6}, {0,10, 4} 


4 









2 





— 


4tt/M 


{0,4+10,0}, {0,6 + 8,0}, 


2 












2 


— 


4tt/M 


{0,0, 4 + 10}, {0,0, 6 + 8} 


2 






2 








Stt/M 


— 


{S + 9,0,0} 


1 


10 


q' 


1 








Att/M 


— 


{16,0,0} 


1 









1 


1 


— 


4tt/M 


{0,4, 12}, {0,6, 10}, {0,8, 8}, 
{0,10, 6}, {0,12, 4} 


S 









2 





— 


4tt/M 


{0,4 + 12,0}, {0,6 + 10,0}, 


2 












2 


— 


4tt/M 


{0,0, 4 + 12}, {0,0, 6 +10} 


2 






2 








Stt/M 


— 


{S + 11, 0,0}, {7 + 9,0,0} 


2 






1 


1 


1 


6tt/M 


Stt/M 


{6,S,S} 


1 


13 



27 



Table 6: The terms through order q^ in the sector m ^ = 'm++ with the macroscopic 

momentum AP = constructed from the rules of sec, 4.3. The minimum momenta 



jDns 

min 



M 



.rrir, 



+ m2s + m_2s + 1) and P, 



±2s 
min 



[m 



+ 



m.23+m_2s+l 



(4.13). Here m2s + rn 



_2s is odd and only m2s > itl_2s are explicitly shown. 



1) are obtained from 



order 


^L 


^2s 


^^~2s 


pns 

min 


p±'2s 
^ min 


|pns^p2s^p-2s| (^^j^i^g Qf ^) 


shift 


states 


total 


q' 





1 





— 


27r/M 


{0,2,0} 





1 


2 


q' 





1 





— 


27r/M 


{0,4,0} 





1 


2 


q' 




1 


1 
1 






Stt/M 


27r/M 
37r/M 


{0,6,0} 
{3,3,0} 






1 
1 


4 


q' 




1 


1 
1 






37i/M 


27r/M 
37r/M 


{0,8,0} 

{3, 5,0}, {5, 3,0} 






1 
2 


6 


q' 




1 


1 
1 






Sn/M 


27r/M 
37r/M 


{0,10,0} 

{3, 7,0}, {5, 5,0}, {7, 3,0} 






1 
3 


8 


<f 




1 




1 

1 

2 





1 


Sn/M 


27r/M 
37r/M 

37r/M 


{0,12,0} 
{3, 9,0}, {5, 7,0} 
{7, 5,0}, {9, 3,0} 
{0,3 + 5,3} 






vr/M 


1 

4 
1 


12 


q' 




1 


2 


1 
1 

2 
1 





1 




37i/M 
A-n/M 


27r/M 
37r/M 

37r/M 
A-k/M 


{0,14,0} 

{3, 11,0}, {5, 9,0} 

{7, 7,0}, {9, 5,0}, {11, 3,0} 

{0,3 + 7, 3}, {0,3 + 5, 5} 

{4 + 6,4,0} 






tt/M 



1 

5 
2 
1 


18 


q' 




1 



2 



1 
1 

2 

1 
3 





1 





37r/M 
Ati/M 


27r/M 
37r/M 

37r/M 

47r/M 
37r/M 


{0,16,0} 

{3, 13,0}, {5, 11,0}, {7, 9,0} 

{9, 7,0}, {11, 5,0}, {13, 3,0} 

{0,3 + 9, 3}, {0,5 + 7, 3} 

{0,3 + 7, 5}, {0,3 + 5, 7} 

{4 + 8, 4,0}, {4 + 6, 6,0} 

{0,3+5+7,0} 






tt/M 



tt/M 


1 

6 

4 
2 
1 


28 



28 



Table 7: The terms through order q^ in the sector m y. = 77i++ with the macroscopic 

momentum AP = n constructed from the rules of sec, 4.3. The minimum momenta 



pns 
min 



j^yiTT'ns 



y±2s 
min 



mV^t. 



+ 



m2s+m-2s+l 



+ m2s + m_2s + 3) and P^,„ - Mv^n^ . 2 

(4.14). Here m2s + 'fn-2s is odd and only m2s > iti^2s are explicitly shown. 



+ 3) are obtained from 



order 




-L 


-1.. 


pns 

min 


p±'2s 
min 


{P"^ P2.^ p-2s| (^j^i^g Qf ^) 


shift 


states 


total 


g' 










— 


4tt/M 


{0,4,0} 





1 


2 


q' 










— 


4tt/M 


{0,6,0} 





1 


2 


q' 










— 


4tt/M 


{0,8,0} 





1 


2 


q' 




1 








Stt/M 


4tt/M 
Stt/M 


{0,10,0} 
{S,S,0} 






1 
1 


4 


q' 




1 








Stt/M 


47r/M 
Stt/M 


{0,12,0} 
{S,7,0},{7,S,0} 






1 
2 


6 


q' 




1 








Stt/M 


4tt/M 
Stt/M 


{0,14,0} 
{S,9,0},{7,7,0},{9,S,0} 






1 
3 


8 


g« 




1 




2 





1 


Stt/M 


4tt/M 
Stt/M 

Stt/M 


{0,16,0} 

{S, 11,0}, {7, 9,0} 

{9,7,0},{11,S,0} 

{0,S + 7,S} 






-tt/M 


1 

4 
1 


12 



Table 8: The terms through order q^ in the sum of m__ 
m_+ = 'm_|_+. These are compared with the terms in q~ 



1/3^1/24^2 ^f (2.30d). 



1 and the AP = term of 



order 


m_+ — m++ = 1 


m_+ — m++ = 


total 


g"V3gl/2452 


q' 


1 





1 


1 


q' 





2 


2 


2 


q' 


1 


2 


3 


3 


q' 


1 


4 


S 


5 


q' 


2 


6 


8 


8 


q' 


S 


8 


13 


13 


q' 


7 


12 


19 


19 


q' 


10 


18 


28 


28 


q' 


13 


28 


41 


41 



29 



Table 9: The terms through order q^ of the sum of m^ 



UIa 



-1 and the AP = tt term 



of m_ 



rrij 



.. These are compared with the terms of q ^/^^g^/^'^62 of (2.30e) 



order 


m_+ — m_|_+ = —1 


m_+ — m++ = 


total 


^-1/12^1/24^2 


g° 


1 





1 


1 


q' 


1 





1 


1 


q' 


1 


2 


3 


3 


q' 


2 


2 


4 


4 


q' 


6 


2 


8 


8 


q' 


7 


4 


11 


11 


q' 


12 


6 


18 


18 


q' 


17 


8 


25 


25 


q' 


26 


12 


38 


38 
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